Adaptive digital filter

ABSTRACT

An adaptive digital filter which is used, for instance, for a prediction filter in an ADPCM modulator and an ADPCM demodulator, for providing stable operation and minimum phase shift response has been found. The filter (FIG. 4, FIG. 5) has three parallel branches (A,B,C), two (A,B) of them have a plurality of series connected non-recursive filter elements (A 1  -A n , B 1  -B n ) each of which has degree not larger than three relating to an operator (Z -1 ), and third branch (C) is merely a conductive line. Tap coefficients (c 1 , c 2 , . . . , d 1 , d 2 , . . . ) of non-recursive filter elements change according to solutions of the transfer function of the filter. A numerator and/or a denominator of the transfer function of the present filter is a Chebychev polynominal relating to an operator (Z -1 ), and has a zero point and/or a pole. The solutions (w i , v i ) of a numerator and a denominator locate alternately on a unit circle (FIG. 6) on a Z -1  plane. When solutions of the transfer function do not locate alternately on said unit circle, tap coefficients (c 1 , c 2 , . . . , d 1 , d 2 , . . . ) are not updated, since non-alternate solutions do not assure stable operation and/or minimum phase shift response of the filter.

BACKGROUND OF THE INVENTION

The present invention relates to an improved iterative controlledadaptive digital filter which is stable in operation and has improvedphase response characteristics. The present filter is used for instanceas a predictive filter in an adaptive PCM (ADPCM) modulator and/ordemodulator. According to the present invention, only the relativelocation of a pair of solutions of a transfer function on a unit circleon a Z⁻¹ plane is monitored, while a prior art monitors the valuesthemselves of the solutions (zero, and pole) of denominator andnumerator of transfer function of a non-recursive digital filter whichhas a feedback loop.

FIG. 1 shows a block diagram of an ADPCM system which is one of theapplicatoins of the present invention. In FIG. 1, the numeral 0 is anA/D converter to convert an analog signal (for instance, voice signal orpicture signal) at the input terminal 11 to a digital form, 1,4 and 7are adders, 2 is a quantizer for converting an input digital signal to aPCM code, 3 and 6 are inverse-quantizers which demodulates a PCM signal,5 and 8 are filters, 9 is a D/A converter for converting a digitalsignal to an analog form, and 11 through 18 are terminals. An analogsignal at the input terminal 11 is converted to a digital form by theA/D converter 0. A digital signal s_(k) (k shows time) is applied to theadder 1 which adds the expected signal s_(k) to the input signal s_(k)and provides the sum which is the error signal e_(k). The error signale_(k) is quantized by the quantizer 2 and is transmitted to an externalcircuit through the output terminal 13. The signal at the terminal 13 isreproduced to an error signal e_(k) by the inverse-quantizer 3 and isapplied to the adder 4 which adds said predicted value s_(k) and the sumis the reproduced value s_(k).

Similarly, the signal at the output terminal 13 is transmitted to areception side or a demodulation side. Noise might be added to thesignal during the signal is transmitted to the reception side. In areception side, the reception signal at the input terminal 16 isreproduced to the reproduce signal s_(k) by the inverse-quantizer 6, theadder 7 and the filter 8. Further, the output signal at the terminal 17is converted to an analog form by the D/A converter 9.

FIG. 2 is a block diagram of a prior non-recursive digital filter forthe filters 5 and 8 in FIG. 1. In FIG. 2, the symbol Z⁻¹ is a delaycircuit which provides the delay time T which is equal to the samplingperiod of the digital signal. The symbols α₁ through α_(n) are tapcoefficients, and the value of them are iteratively adjusted accordingto an input signal so that the error signal e_(k) becomes minimum.

It has been known that an all pole type filter is preferable for aspeech signal. On the other hand, a filter which has not only a pole butalso a zero point is preferable for multi-level digital signal which hasquick change in both amplitude and phase.

However, a filter with both a pole and a zero point has not been usedbecause that kind of filter is apt to oscillate and be affected bynoise.

FIG. 3 shows a block diagram of an ADPCM system which includes both apole and a zero point. In the figure, the numerals 19, 20, 21 and 22 arefilters, 23 through 26 are terminals. Other symbols in FIG. 3 are thesame as those of FIG. 1. The filters 19 through 22 are non-recursivefilters with the structure of FIG. 2.

The transfer function H(Z⁻¹) of the transmission side of FIG. 3 has thefollowing form.

    H(Z.sup.-1)=((1-H.sub.p (Z.sup.-1))/(1+H.sub.z (Z.sup.-1))

     =h.sub.p (Z.sup.-1)/h.sub.z (Z.sup.-1)                    (1)

The transfer function of the reception side is the inverse number of theequation (1). In FIG. 3, the transfer function of the filter between theterminals 12 and 13 is provided a zero point by the non-recursive filter19, and a pole by the non-recursive filter 20. Therefore, the filter 19is called a zero filter, and the filter 20 is called a pole filter. Thereception side has a zero filter 21, and a pole filter 22.

Conventionally, it has been known that a prior filter of FIG. 2 is notstable in operation. Further, in a prior art, the response of the systemis not stable for an input signal which has unexpected statisticalnature, and/or noise, and further, the filter is apt to oscillate and/orthe reproduced code has much code error.

In order to solve the above problem, one solution is to use a fixedfilters 20 and 22 (only filters 19 and 21 are adaptive filters), and theother prior solution is to delete the filters 20 and 22, and that thefilters 19 and 21 have the series connection of a plurality of dualquadratic element filters each of which has a pair of tap coefficientsin a stable area. However, the former solution has the disadvantage thatthe adaptive capability and/or the redundancy compression capability isreduced. The latter solution has the advantage that the stable conditionof the filter is satisfied, but, no mathematical algorithm fordetermining mutual relations between tap coefficients and solutions(pole and zero point) of each element filter. If the filters areadjusted so that an error signal becomes minimum, the convergence ofsolutions becomes slow and effect of a zero point becomes vague.

SUMMARY OF THE INVENTION

It is an object of the present invention to overcome the disadvantagesand limitations of a prior adaptive filter by providing a new andimproved adaptive filter.

It is also an object of the present invention to provide an adaptivefilter which has both a pole and a zero point, stable in operation, andhas excellent phase response.

The above and other objects are attained by an adaptive digital filterwith a transfer function which includes at least one of a pole and azero point in a denominator and a numerator in a transfer function whichchanges successively or at least in every predetermined period,comprising; at least one of denominator and numerator of a transferfunction being a Chebyshev polynominal relating to an operator Z⁻¹ ;said filter having three substantially parallel branches (A,B,C); firstbranch (A) and second branch (B) having a plurality of series connectednon-recursive filter elements which have degree not larger than threerelating to an operator Z⁻¹ ; third branch (C) being a merely conductiveline; and means for determining tap coefficients of said non-recursivefilter elements according to solutions of said transfer function.

Preferably, solutions of a numerator, and solutions of a denominator ofthe transfer function locate alternately on a unit circle on a Z⁻¹ planeto assure stable operation and/or minimum phase shift response of thefilter.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing and other objects, features, and attendant advantages ofthe present invention will be appreciated as the same become betterunderstood by means of the following description and accompanyingdrawings wherein;

FIG. 1 is a prior ADPCM system,

FIG. 2 is a block diagram of a prior adaptive filter,

FIG. 3 is an ADPCM system which the present invention is utilized,

FIG. 4 is a block diagram of the adaptive filter according to thepresent invention when the value (n) is an odd number,

FIG. 5 is a block diagram of the adaptive filter according to thepresent invention when the value (n) is an even number,

FIG. 6 is a graph showing the location of solutions of a transferfunction on a Z⁻¹ plane,

FIG. 7 is a block diagram of another embodiment of the present adaptivefilter when the number (n) is even number,

FIG. 8 is a block diagram of still another embodiment of the presentinvention when the value M is 4,

FIG. 9 is a block diagram of still another embodiment of the presentinvention,

FIG. 10 shows the pattern of the solutions of the transfer function inthe embodiment of FIG. 9,

FIG. 11 shows another pattern of the solutions of the transfer functionin the embodiment of FIG. 9,

FIG. 12 is a block diagram of still another embodiment according to thepresent invention,

FIG. 13 is a block diagram of still another embodiment according to thepresent invention.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

FIG. 4 shows a block diagram of the adaptive filter according to thepresent invention. FIG. 4 is the detailed block diagram of the filter19, 20, 21 or 22 in FIG. 3, and is used by inserting the same in theblock 19, 20, 21 or 22 in FIG. 3. The present filter has threesubstantially parallel branches A, B and C. The former two branches Aand B have a series circuit which has a plurality of non-recursivefilters each of which has the degree equal to or less than 2 relating toor as a function of the operator Z⁻¹. The third branch C is just adirect line. The adder 28 is provided to add the outputs of the branchesA and B, and another adder 29 is provided to add the output of the adder28 and the third branch C. The attenuator 30 for multiplying 1/2 to theinput signal is provided between the input terminal and the branches Aand B. The numeral 27 is another attenuator for providing themultiplication by the value which is 1 or less than 1 but close to 1, tothe input signal. The structure of the filters 19 through 22 arecompletely identical, except that a degree of each element filter, andthe solution (pole or zero point) of the filter depending upon the tapcoefficients are different from those of other filters.

When the filter is a zero filter, each tap coefficient is given below,where n is a degree of an element filter.

    c.sub.1z, c.sub.2z, ,,,, c.sub.nz/2

    d.sub.1z, d.sub.2z, ,,, d.sub.nz/2

When the filter is a pole filter, each tap coefficient is shown below,where m is a degree of an element filter.

    c.sub.1p, c.sub.2p, ,,,, c.sub.mp/2

    d.sub.1p, d.sub.2p, ,,, d.sub.mp/2

Further, it is assumed that each tap coefficient is different from othertap coefficients.

The transfer function h(Z⁻¹) of the present filter is given by thefollowing equations.

When n is an odd number; ##EQU1## When n is an even number; ##EQU2##FIG. 4 shows the case when n is an odd number, and the case when n is aneven number (equation 3) is shown in FIG. 5.

It is assumed for the sake of the simplicity of the explanation that nis an even number, a filter is a zero filter in the followingdescription.

FIG. 6 shows the location of the solutions on a Z⁻¹ plane when n=8. Allthe solutions are imaginary numbers, which are positioned on a unitcircle (|Z⁻¹ |=1). The symbol w_(i) is the solution of a denominator (ora numerator) of a transfer function, and v_(i) is the solution of thenumerator (or denominator) of a transfer function. It should be noted inFIG. 6 that a set of solutions w_(i) and another set of solutions v_(i)are located alternately, that is to say, each solution in a first group(w_(i)) separates a pair of adjacent solutions of another group (v_(i)).Apparently, the sequence of the solutions in FIG. 6 is w₁, v₁, w₂, v₂,w₃, v₃, and therefore, w_(i) and v_(i) appear alternately.

The tap coefficients c_(i), and d_(i) are shown below.

    c.sub.1 =e.sup.jw i=cos w.sub.i

    d.sub.i =e.sup.jv i=cos v.sub.i                            (4)

where w_(i) and v_(i) are angular frequency (radian) of a solution on aunit circle on a Z⁻¹ plane.

The relations between coefficients of a transversal filter in FIG. 2,and coefficients of the present adaptive filter (FIG. 4) are shownbelow, where the transfer function of the transversal filter in FIG. 2is; ##EQU3##

(1) When M=2;

    α.sub.1 =-1/2(c.sub.1 +d.sub.1)                      (6)

    α.sub.2 =1/2(2+c.sub.1 -d.sub.1)                     (7)

(2) When M=4;

    α.sub.1 =-1/2(c.sub.1 +c.sub.2 +d.sub.1 +d.sub.2)    (8)

    α.sub.2 =1/2[4+c.sub.1 +c.sub.2 +c.sub.1 c.sub.2 -(d.sub.1 +d.sub.2 -d.sub.1 d.sub.2)]                                        (9)

    α.sub.3 =-1/2[c.sub.1 +c.sub.2 +c.sub.1 c.sub.2 +d.sub.1 +d.sub.2 -d.sub.1 d.sub.2 ]                                        (10)

    α.sub.4 =1/2[2+c.sub.1 +c.sub.2 -(d.sub.1 +d.sub.2)] (11)

(3) When M=6;

    α.sub.1 =-1/2(c.sub.1 +c.sub.2 +c.sub.3 +d.sub.1 +d.sub.2 +d.sub.3) (12)

    α.sub.2 =1/2[6+c.sub.1 +c.sub.2 +c.sub.3 +c.sub.1 c.sub.2 +c.sub.1 c.sub.3 +c.sub.2 c.sub.3

     -(d.sub.1 +d.sub.2 +d.sub.3)+d.sub.1 d.sub.2 +d.sub.1 d.sub.3 +d.sub.2 d.sub.3 ]                                                 (13)

    α.sub.3 =-1/2[2(c.sub.1 +c.sub.2 +c.sub.3)+c.sub.1 c.sub.2 +c.sub.1 c.sub.3 +c.sub.2 c.sub.3 +c.sub.1 c.sub.2 c.sub.3

     +2(d.sub.1 +d.sub.2 +d.sub.3)-(d.sub.1 d.sub.2 +d.sub.1 d.sub.3 +d.sub.2 d.sub.3)+d.sub.1 d.sub.2 d.sub.3 ]                        (14)

    α.sub.4 =1/2[6+2(c.sub.1 +c.sub.2 +c.sub.3)+c.sub.1 c.sub.2 +c.sub.1 c.sub.3 +c.sub.2 c.sub.3 +c.sub.1 c.sub.2 c.sub.3

     -2(d.sub.1 +d.sub.2 +d.sub.3)+d.sub.1 d.sub.3 +d.sub.1 d.sub.3 +d.sub.2 d.sub.3 -d.sub.1 d.sub.2 d.sub.3 ]                        (15)

    α.sub.5 =-1/2[c.sub.1 +c.sub.2 +c.sub.3 +c.sub.1 c.sub.3 +c.sub.1 c.sub.2 +c.sub.2 c.sub.3

     +d.sub.1 +d.sub.2 +d.sub.3 -(d.sub.1 d.sub.2 +d.sub.1 d.sub.3 +d.sub.2 d.sub.3)]                                                 (16)

    α.sub.6 =1/2[2+c.sub.1 +c.sub.2 +c.sub.3 -(d.sub.1 +d.sub.2 +d.sub.3)](17)

(4) When M=8;

    α.sub.1 =-1/2[c.sub.1 +c.sub.2 +c.sub.3 +c.sub.4 +d.sub.1 +d.sub.2 +d.sub.3 +d.sub.4 ]                                       (18) ##EQU4##

    α.sub.7 =-1/2[c.sub.1 +c.sub.2 +c.sub.3 +c.sub.4 +c.sub.1 c.sub.2 +c.sub.1 c.sub.3 +c.sub.1 c.sub.4 +c.sub.2 c.sub.3 +c.sub.2 c.sub.4 +c.sub.3 c.sub.4

     -(d.sub.1 d.sub.2 +d.sub.3 +d.sub.4 +d.sub.1 d.sub.2 +d.sub.1 d.sub.3 +d.sub.1 d.sub.4 +d.sub.2 d.sub.3 +d.sub.2 d.sub.4 +d.sub.3 d.sub.4)](24)

    α.sub.8 =1/2(c.sub.1 +c.sub.2 +c.sub.3 +c.sub.4 -d.sub.1 -d.sub.2 -d.sub.3 -d.sub.4 +2)                                     (25)

Although in the prior filter of FIG. 2 it is impossible to monitor thestability of the solutions, the present filter has the advantages asfollows since the solutions are located on a unit circle of a Z⁻¹ planeand each solution of a first group separates a pair of adjacentsolutions in another group, and the coefficients of the filter can bedetermined by said equation (4).

(a) The stable operation of a pole filter with any degree is alwaysexpected.

(b) The necessary accuracy of the calculation for the update of asolution and/or a tap coefficient can be reduced as compared with thatof a prior art. In a prior art of FIG. 2, the amplitude and the phase ofthe solutions are directly calculated, and the fact that the solution islocated outside of the circle of FIG. 6 is acknowledged. On the otherhand, in the present invention, the amplitude of all the solutions is 1,and therefore, it is enough to check only a phase of the solutions.

(c) When a solution of a zero filter is simultaneously checked so thatsaid solution is controlled to separate another group of solutions, boththe poles and the zero points of the transfer function are located inthe left half plane of the S-area (S is a Laplace operator), andtherefore, the minimum phase shift response of the circuit is expected.

The successive (iterative) update of the coefficients are accomplishedby for instance correlation calculation, an inverse filtering, etc. Thetransfer function H(Z) of the present successive adaptive filter may beall-pole type, all-zero type, or pole-zero type.

FIG. 7 shows the second embodiment of the present adaptive filter, whenthe degree of the filter is an even number. In the figure, the numeral7-1 is an adaptive filter, 7-2 is an inverse filter for generatinggradient vector component of a filter coefficient d_(i) (i=1-M/2), and7-3 is an inverse filter for generating a gradient vector component of afilter coefficient c_(i) (i=1-M/2). The adaptive filter portion 7-1 inFIG. 7 has at least one branch which has at least one sampling period ofdelay Z⁻¹ between the input 7-4 and the output 7-5, instead of a directline C of FIG. 4. The effect of 7-1 of FIG. 7 is the same as that ofFIG. 4. Each box in FIG. 7 shows Z⁻¹.

According to the algorithm of a gradient method, a filter coefficient ofan adaptive filter at time k+1 is updated by the following equations.

    c.sub.i.sup.k+1 =c.sub.i.sup.k +Δe.sub.k (∂s.sub.k /∂c.sub.i.sup.k)(i=1 - - - M/2)              (26)

    d.sub.i.sup.k+1 =d.sub.i.sup.k +Δe.sub.k (∂s.sub.k /∂c.sub.i.sup.k)(i=1 - - - M/2)              (27)

The gradient vector component of the equations (26) and (27) are shownin the following equations, where E(z) is a z conversion of thedifference signal e_(n).

    ∂s.sub.k /∂c.sub.i.sup.k =-1/2πjφz.sup.k-1 E(z)U.sub.1 (z)(z.sup.-2 -c.sub.i.sup.k z.sup.-1 +1).sup.-1 z.sup.-1 dz (28)

    ∂s.sub.k /∂d.sub.i.sup.k =-1/2πjφz.sup.k-1 E(z)U.sub.2 (z)(z.sup.-2 -d.sub.i.sup.k z.sup.-1 +1).sup.-1 z.sup.-1 dz (29)

The following equations are useful for simplifying a hardware structure,instead of the above two equations.

    c.sub.i.sup.k+1 =c.sub.i.sup.k +Δ·sgn (e.sub.k)·sgn (∂s.sub.k /∂c.sub.i.sup.k)      (30)

    d.sub.1.sup.k+1 =d.sub.i.sup.k +Δ·sgn (e.sub.k)·sgn (∂s.sub.k /∂d.sub.i.sup.k)      (31)

As described above, according to the second embodiment, the gradientvector component of successive updated coefficients according to thegradient method algorithm is obtained by an inverse filter. Therefore,the calculation is simple, and the stability of the updated solutionscan be checked merely by the separation check (alternate location check)of the solutions. The embodiment is applicable to a pole-zero typetransfer function, pole type filter, zero type filter.

As a modification of the second embodiment, the gradient vectorcomponent of a filter coefficient is obtained by correlationcalculation, instead of using an inverse filter.

FIG. 8 is the third embodiment, in which the degree of the filter isM=4. In the figure, the numerals 8-1 and 8-2 are input and outputterminals, respectively, 8-3 is an adaptive filter described in theprevious first or second embodiment, 8-4 is an output buffer of anoutput of the adaptive filter 8-3, and 8-5 is a calculation circuit forthe calculation of the gradient vector component of a filtercoefficient.

The gradient vector component in the equations (26) and (27) is obtainedby the equations (5), and (8) through (11) as follows. ##EQU5##

FIG. 9 shows the fourth embodiment of the present invention. In thefigure, 1 and 2 are prediction filter, 3 is a line, 4 is a memory, 5 and6 are calculation elements. The circuit of FIG. 9 provides theprediction signal s from the input signal s. Then, the difference ebetween the signal s and the prediction signal s is transmitted to theline. Then, the reception side reproduces the reproduction signal s ofthe original signal s. The signal in the reception side is indicated by('), which means that the signal is distorted by noise in thetransmission line. The symbols in the figure are as follows.

s; input pulse signal,

s₁, s₁ '; output of the filter 1

s₂, s₂ '; output of the filter 2

e, e'; difference signal,

s; output signal pulse of the system, that is to say, the reproductionsignal of the signal s.

H₀ ; transfer function of the reception side, input e', output s₀ '

c_(i).sup.(2), d_(i).sup.(2) ; tap coefficients (FIG. 9) of theprediction filter 2.

While the second and the third embodiments control directly the tapcoefficients of an adaptive predictive filter, the fourth embodimentcontrolls the angular frequency w_(i) and v_(i) in the equation (4).

A pair of solutions w_(i).sup.(1) and v_(i).sup.(1) of the predictionfilter 1 in FIG. 9 (i=1 - - - m/2) are controlled according to the nextequation, where k is a discrete time. ##EQU6## The tap coefficients areshown below. ##EQU7## where; ##EQU8## The coefficients A_(i),j.sup.(1)and B_(i),j.sup.(1) when m=6 in the equation (38) are shown below.##EQU9## Similarly, a pair of solutions w_(i).sup.(2), andv_(i).sup.(2), i=1 - - - n/2 are defined by the following equations.##EQU10## The tap coefficients are given below. ##EQU11## where##EQU12##

The constant (α) in the equations (40) and (49) are in the range 2⁻⁷ and2⁻¹⁰.

The four kinds of solutions w_(i).sup.(1), v_(i).sup.(1) (i=1 - - -m/2(prediction filter 1), w_(i).sup.(2), v_(i).sup.(2) (i=1 - - -n/2(prediction filter 2)) are obtained iteratively by the calculation ofthe equations (38) through (49). Those values correspond to the angularfrequency (radian) of the zero point (w_(i).sup.(1), v_(i).sup.(1)), andthe pole (w_(i).sup.(2), v_(i).sup.(2)) of the transfer function H₀.

The fourth embodiment is based upon the fact that the result of theiterative change of the four normalized angular frequency is in apredetermined range, and has the particular pattern defined by a kind ofa signal. FIG. 10 shows the features of a speech, and a modem signal(9.6 kbps). In case of a speech signal, the range of the solutionoverlaps with that of an adjacent solution. On the other hand, in caseof a modem signal, the range of the solutions does not overlap. Thememory 4 stores the upper limit value and the lower limit value of therange of the solution as a window set. The updated values of thesolutions are calculated by the calculation element 5 (equations (38)and (47)). Next, sequence check of the solutions is accomplished. Ifthose conditions are satisfied, the transfer function H₀ (reception),and the inverse transfer function 1/H₀ (transmission side) are stableand provide a minimum phase shift response.

Next, the check if the solutions are within the window set (windowcheck) is accomplished. When the solutions are located in the window setof a modem signal, the input signal is recognized as a modem signal, andthe identification signal ID is output. The sequence check and thewindow check are accomplished by the calculation element 6, which isenough merely to perform subtraction and recognition of a sign (positiveor negative). If the calculation element 6 provides an output signalafter a plurality of signal recognitions, a time constant is obtained ina recognition. When no identification output ID is provided, it isrecognized that a reception signal is not normal, and perhaps a noisedisturbed the signal. In that case, the value of the solution is resetto the previous value or the initial value so that the system does notoscillate. Thus, the unexpected operation by noise is prevented.

As described above, the fourth embodiment transmits a difference of adigital signal, and monitors iteratively the solution of a predictionfilter, then, the stable transmission of a signal is obtained, theminimum phase shift response is satisfied, the recognition of a signalis accomplished both in a transmission side and a reception side, andfurther, the affection by a noise is prevented.

FIG. 11 shows the location of the solutions for a modem signal of 9.6kbps, and a modem signal of 4.8 kbps. In this case, the ranges of thesolutions are similar to each other.

The embodiments of FIG. 10 and FIG. 11 function to discriminate an inputsignal from a speech signal and a modem signal.

FIG. 12 shows a fifth embodiment of the present invention. In thefigure, the symbol AQ is a quantizer, AQ⁻¹ is an inverse quantizer. Thequantizer AQ provides the pulse code train I_(k) from the differencesignal e_(k), and the inverse quantizer AQ⁻¹ functions to inverse saidsignals. In the embodiment, the code train I_(k) is transmitted to aline, and a reproduction signal S_(k) is obtained in the reception side.The symbol T in the figure shows an analog-digital conversion, or adigital-digital speed conversion, and T⁻¹ functions the reverse of T.

According to the embodiment of FIG. 12, the PCM signal which is thedifferential signal is received, and according to said process, thesignal is discriminated correctly in a short time.

FIG. 13 is a block diagram of sixth embodiment of the present invention.The embodiment of FIG. 13 separates the error signal to the error signale_(k) ^(p) by a pole filter, and another error signal e_(k) ^(z) by azero point filter, while the previous embodiments use an error signalwhich is the difference between an input signal and the sum of a firstprediction signal by a pole filter and a second prediction signal by azero point filter, as an objective function for minimizing the error inan iterative control of an adaptive prediction filter. According to theembodiment of FIG. 13, the updated tap coefficient of a pole filter isobtained by the correlation value of e_(k) ^(p) and reproduced outputsignal S_(k), and the updated value of the tap coefficients of a zeropoint filter is obtained by the correlation of e_(k) ^(z) themselves.

The structure of a pole filter and a zero point filter in FIG. 13 is thesame as that of the embodiments of FIGS. 1 through 5, and thecalculation process in FIG. 13 is shown below.

The control is accomplished according to the following equation at adiscrete time k.

(a) a zero point filter; ##EQU13## The tap coefficients are obtained bythe following equations. ##EQU14## where; ##EQU15## The value (δ_(i)) isusually in the range between 2⁻⁶ and 2⁻¹¹, and the value (α) is in therange between 2⁻⁷ and 2⁻¹⁰. The values A_(i),j.sup.(1) andB_(i),j.sup.(1) in the equation (50) are, when m=6, the same as those ofequations (41) through (46) in the fourth embodiment.

(b) a pole filter;

Similar to the case of a zero point filter; ##EQU16## The tapcoefficients are; ##EQU17## where; ##EQU18##

From the foregoing, it will now be apparent that a new and improvedadaptive filter has been found. It should be understood of course thatthe embodiments diclosed are merely illustrative and are not intended tolimit the scope of the invention. Reference should be made to theappended claims, therefore, rather than the specification as indicatingthe scope of the invention.

What is claimed is:
 1. An adaptive digital filter with a transferfunction which includes at least one of a pole and a zero point in adenominator and a numerator in a transfer function which changessuccessively or at least in every predetermined period, comprising;atleast one of the denominator and numerator of a transfer function in thecomplex Z plane being a Chebyshev polynominal, said filter having threebranches (A,B,C), first branch (A) and second branch (B) having aplurality of series connected non-recursive filter elements (A₁ -A_(n),B₁ -B_(n)) in which said filter elements include delay circuits, a tapand an adder and have degree not larger than three, third branch (C)being a merely conductive line of said second branch of non-recursivefilter elements, branch (C) being coupled with an input terminal of thefilter to which said first branch of branches (A) and (B) are coupledthrough an attenuator with a multiplier 1/2, outputs of those threebranches being coupled together with an output terminal of the filter,and means for determining tap coefficients (c₁, c₂, (B₁, B₂, . . . ),and tap coefficients (d₁, d₂, . . . ) of said non-recursive filterelements (A₁, A₂, . . . ) according to each coefficient of saidChebyshev polynominal.
 2. An adaptive digital filter according to claim1, wherein an inverse filter is provided for providing a tap coefficientof said non-recursive filter elements.
 3. An adaptive digital filteraccording to claim 2, wherein said means for determining tapcoefficients comprises;means for providing an updated tap coefficient bysubtracting at least a part of gradient signal obtained by an inversefilter of said non-recursive filter elements, from previous tapcoefficient, means for monitoring location of solutions on a Z⁻¹ planeif updated solutions for a zero point, and updated solutions for a polelocate alternately on a unit circle on said Z⁻¹ plane to assure at leastone of stable operation of the filter and minimum phase shift responseof the filter, means for re-updating tap coefficient when said solutionsfor a zero point and a pole do not locate alternately on said unitcircle so that those solutions locate alternately, to determine tapcoefficients relating to re-updated solutions.
 4. An adaptive digitalfilter according to claim 1, wherein a check circuit for monitoringrange of updated solutions, and recognizing a kind of an input signal tothe present adaptive digital filter.
 5. An adaptive digital filteraccording to claim 1, wherein means for limiting noise effect isprovided, and said means monitors updated solutions, and when an updatedsolution is out of predetermined range, said means resets a solution toan initial value, or return to a previous value of solution.